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ABSTRACT 


The  problem  of  extracting  the  Rayleigh  phase  of  an  event  of 
interest  from  a  time -overlapping  Rayleigh  phase  of  another  event  was  studied. 
Various  configurations  of  the  long-period  LASA  sensors  and  the  following  pro¬ 
cessing  schemes  were  investigated: 

•  Beamsteer  processing  using  the  great  circle  route 
and  the  dispersion  curve  velocity  occurring  at  the 
peak  power  frequency 

•  Wiener  signal-extraction  multichannel- filter 
processing 

•  A  nonconventional  adaptive  filter  processing 
scheme 

The  following  conclusions  are  based  on  the  results  of  this  work: 

•  Current  frequency -wavenumber  spectral  techniques 
are  not  useful  in  providing  location  information 
needed  for  processor  design 

•  Despite  the  appreciable  differences  in  waveform  as 
seen  across  the  array,  neither  the  beamsteer  nor 
MCF  processors  introduce  serious  distortion  of 
the  target  event 

•  As  expected,  suppression  of  the  interfering  event 
becomes  less  as  the  azimuthal  separation  between 
events  is  reduced  (best  suppression  achieved: 

A=55°,  18  db;  A=28.2°,  12  db) 

•  The  MCF  processor  is  as  much  as  12  db  superior  to 
the  beamsteer  processor  .n  suppressing  the  inter¬ 
fering  event 

i 

•  Multichannel  filtering  of  small-array  configurations 
such  as  the  AO  and  D  ring  or  AO  and  C  ring  is  the 
most  effective  processing  method  of  those  considered 

•  The  nonconventional  adaptive  processor  is  not 
effective  in  solving  this  problem 
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SECTION  I 
INTRODUCTION 

Various  methods  were  studied  for  separating  Rayleigh  waves 
from  two  distinct  epicenters  received  simultaneously  by  a  long-period  array. 
The  problem  was  treated  in  an  off-line  context;  thus,  it  is  assumed  that  both 
events  have  been  detected  and  located  by  other  means.  The  problem  is  to 
process  the  received  data  so  that  one  event  is  adequately  suppressed  for  those 
parametere  useful  for  classifying  the  other  event  to  become  available. 

The  usefulness  of  various  types  of  frequency-wavenumber  (f-k) 
spectra  in  accurately  determining  the  apparent  locations  of  the  two  events  in 
wavenumber  space  is  detailed  in  Section  II.  Knowledge  of  the  apparent  event 
locations  is  important  in  cases  where  significant  lateral  refraction  of  the  sur¬ 
face  waves  occurs  and  tho  great  circle  route  becomes  a  poor  estimate  of  the 
propagation  vector.  In  Sections  III  and  IV,  details  of  the  following  three  types 
of  signal  processors  are  presented: 

•  Beamsteer  processing  using  the  great  circle  route 
and  the  dispersion  curve  velocity  occurring  at  the 
peak  power  frequency 

•  Wiener  signal-extraction  multichannel  filter 
processing 

•  A  nonconventional  adaptive  filter  processing 
scheme 

At  the  time  this  work  was  performed,  time -overlapping  events 
were  not  available.  Simulation  of  the  problem  was  accomplished  by  composit¬ 
ing  individual  Rayleigh- wave  arrivals.  Pertinent  details  of  the  events  used 
are  given  in  Table  1-1. 
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SECTION  II 

ESTIMATION  OF  EVENT  VELOCITY  USING  FREQUENCY- 

WAVENUMBER  SPECTRA 

The  various  types  of  MCF  processors  required  for  separating 
overlapping  Rayleigh- wave  events  require  modeling  of  either  one  or  both 
interfering  events.  This  necessitates  knowing  accurately  the  events'  velocity 
and  direction  of  propagation.  Using  the  great  circle  route  from  epicenter  to 
array  for  direction  and  the  best  current  estimates  of  dispersion  curves  in  ihe 
vicinity  of  the  array  for  velocity  is  the  most  straightforward  approach.  How¬ 
ever,  it  is  well  known  that  lateral  refraction  of  Rayleigh  waves  occurs. 
Consequently,  the  straightforward  estimate  of  event  location  in  frequency- 
wavenumber  space  may  be  a  poor  one.  It  will  be  shown  later  in  this  report, 
for  example,  that  the  apparent  propagation  vector  at  LASA  for  one  event 
originating  in  California  differs  by  about  5°  from  the  great  circle  path. 

—4 

The  usefulness  of  various  types  of  f-k  spectra  in  obtaining 
better  location  estimates  was  investigated.  Everts  from  California  and 
Greenland  were  composited  to  simulate  overlapping  events;  approximate 
Greenlan  1-to-California  amplitude  ratios  of  2:1,  5:1,  and  10:1  were  used. 

«-4 

The  following  types  of  f-k  spectra  were  computed  for  each  of  the  three  com¬ 
posite  ratios: 

•  Conventional  -  conventional  f-k  spectra 

•  HR  1  -  high- re  solution  f-k  spectra  with  the  AO 

sensor  used  as  reference 

/ 
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•  HR  2  -  high- resolution  f-k  spectra  with  the  F- 
ring  sensors  used  as  reference 

Analytic  details  of  these  techniques  are  reported  in  Special  Scientific  Report 

1-4 

No.  2.  The  f-k  spectra  were  calculated  at  0.04  Hz  which  is  near  the  peak 
of  the  single-channel  power-density  spectra  for  both  of  these  events.  The 
composite -event  time  gate  used  for  calculating  the  f-k  spectra  was  330  sec 
and  contained  most  of  the  signal  energy  for  both  events. 
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Proper  evaluation  of  the  f-k  spectral  techniques  requires  ac- 
curate  knowledge  of  the  apparent  f-k  space  location  of  the  events  as  they  pro¬ 
pagate  across  the  array.  Table  II- 1  lists  two  different  reference  estimates 
of  these  locations  as  well  as  those  obtained  from  the  various  f-k  spectra. 

The  computer  program  used  in  generating  these  spectra  outputs  in  the  form 
of  relative  db  power  levels  at  various  points  in  wavenumber  spaces  as  defined 
by  a  linear  rectangular  grid.  The  grid  consists  of  59  locations  in  the  y  direc¬ 
tion  and  49  locations  in  the  x  direction,  both  beginning  in  the  lower  left-hand 
corner.  The  locations  in  Table  II- 1  are  given  in  terms  of  these  coordinates. 
The  first  estimate  was  obtained  using  the  great  circle  routes  from  the  PDF) 
epicenters  to  the  array  and  the  phase  velocity  from  a  LASA  dispersion  curve 
at  0.04  Hz  (first  line,  Table  II- 1).  The  second  estimate  was  obtained  from  a 
previous  study  of  these  two  events  aimed  at  experimental  determination  of  the 
LASA  Rayleigh-wave  dispersion  curve.  ^  This  work  provided  estimates  of  the 
velocity  and  direction  for  both  of  these  events  using  various  3- sensor  combina- 
ions  from  the  LASA.  The  apparent  f-k  space  locations  for  the  two  events 
have  been  visually  estimated  from  these  results  and  are  given  in  the  second 
line  of  Table  II- 1.  The  two  estimates  of  the  Greenland  location  are  in  rela¬ 
tively  close  agreement,  but  the  apparent  California  direction  departs  signifi¬ 
cantly  from  the  great  circle  route. 

Examples  of  the  contoured  f-k  spectra  obtained  are  given  in 
Figure  II- l.  These  are  the  three  different  types  of  f-k  spectra  for  the 
composite  event  using  a  2:1  amplitude  ratio.  The  location  of  peak  power  is 
indicated  by  a  plus  (+)  sign,  and  the  contburs  give  power  levels  in  db  relative 
to  this  peak  power. 
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Greenland/California  =  10:1 
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As  expected,  the  location  of  peak  power  on  all  the  f-k  spectra 
is  near  the  larger  of  the  two  events  (Greenland).  Estimates  of  the  Greenland 
location  (listed  in  Table  II-  1)  are  these  peak-power  locations.  Assuming  that 
the  reference  estimates  of  the  Greenland  event  are  correct,  the  conventional 
and  HR2  spectra  pxoperly  locate  this  event  for  ratios  of  5:1  and  10:1.  For  a 
ratio  of  2:1,  the  Greenland  event  is  improperly  located  by  all  three  types  of 

—4 

f-k  spectra.  For  all  three  types  of  spectra,  the  presence  of  a  relatively 
strong  secondary  event  distorts  the  location  of  the  strongest  event. 

In  all  the  spectra  computed,  a  secondary  peak  occurs  in  the 
vicinity  of  the  California  location.  Since  the  computer  program  does  not 
provide  the  coordinate  where  the  largest  value  of  this  peak  occurs,  the 
coordinate  nearest  its  center  was  selected  as  the  California  estimate.  These 
estimates  differ  seriously  from  either  reference  location  and  must  be  con¬ 
sidered  erroneous. 

If  the  problem  is  to  extract  a  small  event  from  a  much  larger 
interfering  event,  the  results  imply  the  following  conclusions. 

•  It  is  probably  possible  to  obtain  a  good  estimate  of 
the  wavenumber  location  for  the  large  event  using 
conventional  or  HR2  spectra 

•  If  severe  lateral  refraction  of  the  event  has  occurred, 
this  estimate  will  be  superior  to  the  great  circle 
estimate 

•  It  does  not  appear  possible  with  any  of  these 
techniques  to  obtain  an  accurate  estimate  of 
the  weak- event  location 

I 

In  Section  III  it  is  seen  that  the  best  MCF  results  are  obtained  by  using  mea¬ 
sured  data  to  characterize  the  large  interfering  event  and  by  modeling  the 
small  event.  Therefore,  accurate  location  of  the  large  event  in  the  f-k  space 
is  of  little  value.  In  view  of  this,  it  appears  that  none  of  these  f-k  spectra 
contribute  significantly  to  solving  the  interfering- event  problem. 
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For  the  work  discussed  in  this  section,  it  is  assumed  that  the 
Rayleigh  wave  from  a  large  natural  event  has  been  received  and  that  the 
Rayleigh  wave  from  an  event  of  interest  has  been  received  simultaneously. 

It  is  assumed  that  both  events  have  been  detected  and  that  their  epicenters 
have  been  located  by  other  means.  The  problem  is  to  suppress  the  interfering 
event  sufficiently  to  determine  classification  parameters  of  the  event  of 
interest. 

This  problem  has  been  studied  using  two  types  of  processors: 
time- shift-and- sum  and  conventional  MCF.  A  third  approach  is  discussed  in 
Section  IV  of  this  report.  Specific  questions  to  be  answered  here  are  as  follows 

•  How  do  the  two  approaches  compare? 

•  How  does  performance  depend  on  the  azimuthal 
separation  of  the  two  events? 

•  What  configuration  of  the  LASA  long -period 
sensors  works  best? 

•  How  much  suppression  of  the  interfering  event 
can  be  achieved  using  these  methods? 

At  the  time  this  processing  began,  no  known  instances  of  over¬ 
lapping  long-period  events  existed  in  the  available  data  library.  To  simulate 
the  problem,  therefore,  it  was  necessary  to  scale  and  add  individual  events. 
Pertinent  details  of  the  events  used  are  given  in  Table  I-i.  Th»  first  com- 

I 

posite  event  studied  consists  of  the  interfering  Hokkaido  event  and  the  target 
New  Hebrides  event  with  an  azimuthal  separation  of  55°.  Before  adding  these 
t  /o  events,  the  New  Hebrides  data  was  scaled  down  to  make  the  peak  power 
level  of  its  AO  seismometer  output  17  db  below  that  for  the  Hokkaido  event. 
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Since  both  the  beamsteer  and  MCF  processors  were  linear  and 
the  "signal"  and  "noise"  events  were  available  individually,  evaluation  of  the 
processors  consisted  of  applying  the  processors  to  the  uncomposited  individual 
events.  This  allowed  a  more  precise  estimate  of  the  processors'  performance. 
Outputs  of  various  beamsteer  processors  for  the  New  Hebrides  event  were 
illustrated  in  Special  Scientific  Report  No.  20.  A  surface-wave  velocity  of 
3.5  km/sec  was  used  in  computing  all  beamsteer  time  delays  in  this  section. 

In  view  of  the  New  Hebrides  event's  variability,  as  recorded  by  the  various 
vertical  sensors,  these  results  show  a  surprisingly  small  degree  of  signal 
degradation.  The  result  of  steering  various  combinations  of  the  Hokkaido 
traces,  using  the  New  Hebrides  azimuth,  is  shown  in  Figure  III  - 1 .  Note  that 
all  New  Hebrides  traces  have  been  scaled  up  by  a  factor  of  five  in  this  figure 
for  clarity  in  the  display.  Despite  the  large  azimuthal  separation  of  these 
two  events,  the  best  beamsteer  processor  provides  only  6  db  suppression  of 
the  peak  value  of  the  interfering  wave. 

Two  approaches  were  considered  when  designing  the  MCF 
processors.  In  one  case,  both  the  signal  and  noise  correlation  matrices  were 
obtained  from  theoretical  dispersive  plane-wave  models  propagating  along  the 
great  circle  pa*hs  to  LASA  from  New  Hebrides  and  Hokkaido,  respectively. 

The  frequency-dependent  propagational  velocities  used  in  the  models  were 
those  of  the  LASA-North  dispersion  curve.  The  single-sensor  power- 
density  spectrum  for  both  models  had  the  shape  given  in  Figure  III-2.  This 
was  chosen  to  force  the  MCF's  to  work  in  the  frequency  range  where  signal 
and  noise  power  are  known  to  exist.  A  signal-to-noise  ratio  of  approximately 
4  was  used  in  all  MCF  designs.  Two  theoretical  MCF  processors  were 
designed  for  the  New  Hebrides -Hokkaido  composite,  one  employing  the  AO 
and  C-ring  verticals  and  the  other  employing  the  AO  and  E-ring  verticals. 
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FREQUENCY  (H«) 

Figure  III-2.  Theoretical  Model  Power-Density  Spectrum 


In  the  second  type  of  MCF  design,  a  measured  correlation  set 
was  obtained  from  the  composite  traces.  The  fact  that  this  set  contains  a 
small  amount  of  signal  data  was  neglected  and  this  set  was  used  as  the  noise 
correlations.  Signal  correlations  were  again  obtained  from  a  theoretical 
model,  but  in  this  case,  the  single-sensor  power-density  spectrum  was  given 
the  same  shape  as  the  AO  sensor-measured  power-density  spectrum.  An 
MCF  of  this  type  was  designed  for  the  AO  and  C-ring  sensors.  In  this  design, 
the  signal  autocorrelations  were  scaled  by  1.  01  representing  statistical  gain 
fluctuation,  and  the  zero  lags  of  the  noise  autocorrelations  were  scaled  by  1.01 
representing  white  noise. 

The  results  of  applying  these  MCF  processors  to  the  individual 
Hokkaido  and  New  Hebrides  data  sets  are  given  in  Figure  III- 1 .  Examining 
the  data  in  Figure  UI-1  indicates  that 

•  Neither  beamsteer  nor  MCF  processing  introduces 
severe  signal  distortion 

i 

•  Both  AO  and  C-ring  MCF  processors  suppress  the 
peak  value  of  the  noise  by  about  18  db  and  are  sig¬ 
nificantly  superior  to  any  of  the  beamsteer  proces¬ 
sors  in  this  respect 

•  If  the  AO/E-ring  combination  is  used,  neither 
beamsteer  nor  theoretical  MCF  processing  is 
clearly  superior 

•  The  theoretical  AO/E-ring  MCF  is  inferior  to 
either  AO/C-ring  MCF  processor 


Maximum-entropy  power-density  spectra  of  the  MCF  outputs 
and  of  the  AO  Hokkaido  and  New  Hebrides  traces  are  illustrated  in  Figure 
IU-3.  The  C-ring  MCF's  provide  significant  suppression  of  the  Hokkaido 
event  throughout  the  frequency  range  where  appreciable  power  exists.  The 
E-ring  processor,  however,  is  largely  ineffective  beyond  0.05  Hz. 

Frequency- wavenumber  responses  of  the  two  theoretical  MCF 
processors  are  shown  in  Figures  III -4  and  III- 5 .  Only  that  portion  of  wave- 
number  space  which  includes  the  apparent  location  of  the  two  events  is  dis¬ 
played.  The  indicated  frequencies  bracket  the  range  of  peak  power  for  these 
events.  At  each  frequency,  the  points  corresponding  to  the  velocity  of  the 
LASA-North  dispersion  curve  and  to  the  directions  of  the  two  events  are 
indicated  by  dots.  In  all  cases>  the  response  is  very  nearly  0  db  at  the  signal 
location  and  shows  varying  amounts  of  suppression  at  the  noise  location.  At 
0.  07  Hz, the  E-ring  MCF  response  has  a  sharp  gradient  in  the  vicinity  of  the 
noise.  The  poor  performance  of  this  MCF  is  probably  either  due  to  one  or  a 
combination  of  two  causes.  In  view  of  the  snarp  response  gradient  at  0.  07  Hz, 
small  errors  in  locating  the  noise  model  could  result  in  poor  suppression  of 
the  actual  noise.  Alternately,  the  wide  variations  in  the  noise  waveform  — 
as  perceived  by  the  E-ring  sensors  —  may  prevent  effective  MCF  processing 
with  this  configuration.  This  approach  was  discarded  for  subsequent  designs , 
since  using  a  theoretical  noise  model  allows  the  possibility  of  "missing"  the 
actual  interfering  event  and  since  comparison  of  the  two  C-ring  designs  shows 

no  apparent  advantage  for  the  theoretical  design. 

! 

A  second  composite  event,  having  an  azimuthal  separation  of 
28.2  ,  was  formed  from  the  Hokkaido  and  Mongolia  events.  Again,  the 
Hokkaido  event  was  considered  to  be  the  interfering  event,  and  the  Mongolia 
data  was  scaled  down  to  make  its  A0  seismometer  peak  power  level  18  db 
below  that  for  the  Hokkaido  event. 
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Figure  III- 3 .  Maximum- Entropy  Power-Density  Spectra 
of  MCF  Outputs 
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Beamsteer  outputs  for  the  Mongolia  event  were  illustrated  in 

•j 

Special  Report  No.  20.  Again,  signal  degradation  resulting  from  beamsteer 
processing  is  seen  to  be  negligible.  Note  that  about  the  first  third  of  these 
beamsteer  traces  represents  energy  arriving  prior  to  the  onset  of  the  Rayleigh 
phase.  The  results  of  using  the  Mongolia  azimuth  to  steer  various  combina¬ 
tions  of  the  Hokkaido  traces  are  shown  in  Figure  III  -  6 .  The  best  beamsteer 
suppresses  the  peak  value  of  the  interfering  wave  by  less  than  4  db. 

Initially,  two  MCF  processors  were  designed  —  an  11 -point 
filter  using  the  AO  and  C-ring  s»  nsors  and  a  31 -point  filter  using  the  AO  and 
E-ring  sensors.  Both  statistical  gain  fluctuation  and  white  noise  were  used 
in  these  designs.  Noise  and  signal  outputs  for  these  processors,  as  well  as 
those  to  be  discussed  below,  appear  in  Figure  III- 6 .  A  small  degree  of  signal 
degradation  occurs  in  each  case.  Attenuation  of  the  interfering  event  by  these 
two  processors  is  similar  and  is  somewhat  better  than  that  provided  by  any  of 
the  beamsteer  processors.  Nevertheless,  the  peak  value  of  either  MCF  noise 
output  is  only  about  9  db  below  that  of  the  unprocessed  interfering  event  as 
seen  by  the  AO  seismometer.  Apparently,  the  C-ring  aperture  is  insufficient 
to  adequately  resolve  these  two  events,  while  the  E-ring  MCF  suffers  from 
dissimilarities  in  the  interfering  event  as  seen  by  its  sensors. 

In  an  effort  to  permit  the  AO/C-Ring  MCF  greater  latitude  in 
working  on  the  interfering  event,  the  design  was  repeated  without  statistical 
gain  fluctuation.  The  results  in  Figure  III -6  show  no  significant  improvement. 

Power  spectra  of  these  three  MCF  outputs,  along  with  those  of 

t 

the  AO  Hokkaido  and  Mongolia  traces,  are  shown  in  Figure  IU-7.  In  all  cases, 
the  peak  noise  suppression  occurs  at  about  0.  042  Hz.  Its  value  is  13  db  for 
the  ^  ring  MCF  designed  with  statistical  gain  fluctuation  and  15  db  for  the 
other  t  vc  designs.  Average  noise  suppression  over  the  0.03"  to0.07-Hz 
range  is  somewhat  less. 
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Figure  III-6.  Mongolia- Hokkaido  Processing  Results 
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Two  additional  MCF  processors  were  designed  for  this  pair 
of  events.  The  first  was  a  9-channel  25-point  MCF  using  the  AO,  C-,  and  D-ring 
sensors.  The  second  employed  a  combination  of  vertical  and  horizontal  sensors. 
Before  processing,  the  N-S  and  E-W  horizontal  traces  were  scaled  and  trigonometri¬ 
cally  rotated  to  form  a  resultant  trace  which  would  be  seen  by  a  sensor  oriented 
transverse  to  the  Mongolian  signal’s  direction  of  propagation.  An  8 -channel 
11-point  MCF  was  designed  using  the  AO  and  C-ring  vertical  sensors  and  the 
A0f  B2,  and  B3  transverse  horizontal  sensors.  (The  B1  and  B4  horizontal  sensors 
were  anomalous  for  this  composite). 

The  9-channel  (A-,  C-,and  D-ring  vertical  sensors)  processor 
provides  only  a  small  amount  of  improvement  in  noise  attenuation  and  has  rather 
severe  signal  degradation.  The  8-channel  (vertical/horizontal  sensor)  processor 
likewise  suppresses  the  interfering  event  slightly  better  than  the  other  processors 
but  shows  very  little  singal  attenuation. 

With  the  exception  of  the  9-channel  AO,  C-,  and  D-ring  MCF, 
all  processors  studied  have  produced  only  a  small  amount  of  signal  degradation. 
Therefore,  the  figure  of  merit  used  in  summarizing  these  results  will  be  the 
ratio  of  the  peak  amplitude  in  the  processor  noise  output  to  the  peak  amplitude 
in  the  AO  sensor  noise  trace.  Using  this  yardstick,  the  following  conclusions 
can  be  made. 

•  MCF  processing  significantly  suppresses  the 
interfering  wave  better  than  beamsteer  processing. 

The  degree  of  this  difference  varies  with  the 
separation  of  events  and  with  the  sensor  configur¬ 
ation  and  ranges  as  high  as  12  db.  The  best  MCF 
results  are  obtained  by'  using  the  received  composite 
waveforms  to  obtain  an  estimate  of  the  noise 
correlation  matrix  for  filter  design. 

•  In  comparing  5-channel  MCF  processors,  the  AO/ 

C-ring  combination  appears  to  work  as  well  or 
better  than  the  AO/E-ring.  Upon  adding  D-ring 
sensors  to  the  AO,  the  C-ring  combination  substantially 
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increases  the  signal  degradation.  Adding  several 
horizontal  transverse  sensors  to  the  AO/C-ring 
combination  adds  slightly  to  the  suppression 
capability,  with  no  noticeable  increase  in  signal 
degradation. 

•  The  best  noise  suppression  achieved  for  the  events 
separated  by  55°  was  about  18  db.  When  the  event 
separation  was  28.2°,  this  suppression  dropped 
to  12  db. 
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SECTION  IV 

SEPARATION  OF  EVENTS  USING  ADAPTIVE  PROCESSING 

In  the  ideal  case,  two  interfering  plane  waves  can  be  separated 
using  only  two  channels  of  data.  In  practice,  however,  departures  from  plane- 
wave  theory  caused  by  such  phenomena  as  local  scattering  and  multipath  pro¬ 
pagation  seriously  hamper  this  simple  approach.  A  straightforward  adaptive 
filtering  scheme  has  been  studied  in  an  effort  to  overcome  these  difficulties. 

Again,  it  is  assumed  that  two  time -overlapping  Rayleigh  waves 
have  been  received  and  that  the  azimuth  from  station  to  epicenter  for  both  is 
known.  In  this  study,  events  from  Greenland  and  California  were  composited, 
with  the  latter  considered  to  be  the  signal.  In  forming  the  composite,  the 
California  event  was  scaled  down  to  make  the  peak  value  of  its  AO  trace  about 
15  db  below  that  of  the  Greenland  event.  The  composite  has  two  significant 
characteristics: 

•  The  similarity  of  these  events  as  seen  by 
the  various  LASA  long-period  vertical  seis¬ 
mometers  is  better  than  that  observed  for 
most  Rayleigh- wave  arrivals 

•  The  azimuthal  separation  of  these  two  events 
is  approximately  125 

As  a  result,  this  particular  composite  represents  an  unusually  simple 
problem,  and  the  results  obtained  should  provide  an  upper  performance  limit 
for  the  method  under  study. 

After  forming  the  composite  event,  a  beam  was  formed  toward 
the  Greenland  epicenter  using  a  3.  5 -km/sec  velocity.  Presumably  this  beam- 
steer  output  is  a  good  replica  of  the  Greenland  arrival  and  contains  little 
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California  energy.  After  time -shifting  the  individual  composite  traces  to  make 
their  Greenland  arrivals  time-align  with  that  in  the  beamsteer,  individual  short 
1 -channel  adaptive  prediction  filters  were  designed  to  predict  each  of  the  traces 
from  the  beamsteer  output.  Ideally,  each  of  the  prediction-error  traces  should 
contain  the  undistorted  California  event  and  the  strongly  attenuated  Greenland 
event.  Finally,  the  prediction-error  traces  were  beamsteered  for  the  Cali¬ 
fornia  azimuth  to  further  suppress  the  remaining  Greenland  energy. 

The  adaptive  algorithm  used  in  designing  the  single-channel 

4 

prediction  filters  was  given  previously  and  is  repeated  here. 


«  =  y  -  F  X 

n  7n  n  n 


F  =  F  +  2k  €  X 
n+ 1  n  s  n  n 

where 

y^  is  the  n1*1  data  point  in  the  trace  to  be  predicted 

is  an  N  x  1  column  vector  containing  the  N  points 
of  the  beamsteer  output  used  in  predicting  y^ 

F^  is  an  N  x  1  column  vector  containing  the  N  filter 
weights  used  in  predicting  y^ 

€  ^  is  the  error  in  predicting  y^ 

k^  is  a  convergence  parameter  determining  the  rate 
of  filter  adaptation 

Widrow  has  shown  that,  for  an  N-point  adaptive  filter  which  predicts  a  time- 
series  data  point  from  the  preceding  N  data  points,  the  following  stability  re 
quirement  holds  for 


0  <  k  < 


T 

X  X 
jmax  jmax 
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In  the  present  case  where  one  time  series  is  used  to  predict  another  similar 
time  series,  a  practically  identical  constraint  is  as  follows: 


0  <  k  <  — - 

x  (NFP) 
jmax 


where 


x  is  the  maximum  value  of  the  data 

jmax 

NFP  is  the  number  of  filter  points 


For  the  type  of  data  used,  this  upper  limit  on  kg  is  conservative.  Since  it  is 

desirable  to  have  the  filters  adapt  as  rapidly  as  possible,  tht  beamsteer  trace 

is  first  scanned  to  determine  x.  ,  and  k  is  set  to  the  upper  limit.  The 

jmax  s 

use  of  short  filters  is  necessary  to  permit  sufficiently  rapid  adaptation. 

The  technique  was  first  evaluated  using  the  AO,  D-,  and  E-ring 
sensors  with  3-point  filters.  Figure  IV  *1  shows  the  AO -sensor  California  and 
Greenland  traces,  the  nine  prediction-error  traces,  and  the  final  output.  Note 
that  the  Greenland  AO  trace  has  been  scaled  down  by  a  factor  of  10  in  this  display, 
Clipping  of  several  of  the  prediction-error  traces  is  a  plotting  error. 

The  first  portion  of  the  interfering  event  is  attenuated  by  more 
than  30  db;  however,  as  the  event  proceeds,  this  attenuation  diminishes  rapidly 
and  only  about  10-db  suppression  of  the  interfering  event  is  observed  near  the 
end  of  the  record.  Apparently,  dissimilarity  in  the  individual  seismometer 
outputs  (which  becomes  more  severe  in  the  latter  portion  of  the  record)  seriously 
hampers  the  single -channel  filters.  In  addition  to  this  limited  effectiveness  in 
suppressing  the  interfering  event,  about  6-db  signal  attenuation  is  observed. 
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Figure  IV-  1.  AO,  D-, and  E-Ring  Adaptive  Processor  Results 
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Figure  IV-2.  AO,  C-,and  D-Ring  Adaptive  Processor  Results 
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In  an  effort  to  lessen  the  effect  of  signal  dissimilarities,  the 
experiment  was  repeated  using  the  AO,  C-,  and  D-ring  sensors.  In  this  case, 

1 -point  filters  were  used,  permitting  a  threefold  increase  in  the  rate  of 
adaptation  under  the  previously  cited  stability  constraint.  The  results  in  this 
case  are  shown  in  Figure  IV-2.  Noise  attenuation  in  the  latter  portion  of  the 
record  has  improved  only  slightly.  Similarly,  degradation  of  the  signal  has 
lessened  but  is  still  severe.  In  addition,  suppression  of  the  first  portion  of 
the  interfering  event  is  reduced  to  approximately  24  db. 

Despite  the  relatively  favorable  nature  of  this  particular 
compound  event,  the  technique  appears  to  be  very  limited  in  its  ability  to 
separate  the  two  events.  Most  cases  of  practical  interest  will  be  less  favorable* 
both  from  the  standpoint  of  azimuthal  separation  and  of  signal  dissimilarity 
across  the  array.  Therefore,  it  is  concluded  that  the  technique  is  of  little 
value  for  separating  time -overlapping  long-period  Rayleigh  waves. 


SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  following  conclusions  are  based  on  the  results  of  this  work. 

•  Current  frequency  wavenumber  spectral 
techniques  are  not  useful  in  providing  location 
information  needed  for  processor  design 

•  Despite  the  appreciable  differences  in  waveform 
as  seen  across  the  array,  neither  the  beamsteer 
nor  MCF  processors  introduce  serious  distortion 
of  the  target  event 

•  As  expected,  suppression  of  the  interfering 
event  becomes  less  as  the  azimuthal  separation 
between  events  is  reduced  (best  suppression 
achieved:  A  =  55°,  18db;A=  28  2°,  12  db) 

•  The  MCF  processor  is  superior  to  'he  beamsteer 
processor  by  as  much  as  12  db  in  suppressing  the 
interfering  event 

•  Multichannel  filtering  of  small-array  configurations 
(e.g.  ,  the  AO/D  ring  or  AO/C  ring)  is  the  most  ef¬ 
fective  processing  method  of  thOv>e  considered 

•  The  nonconventionai  adaptive  processor  is  not 
effective  in  solving  this  problem 

Further  study  is  recommended  in  the  following  areas. 

•  One  example  of  the  use  of  horizontal  sensors  in 
event  separation  has  been  investigated  and 
further  exploration  seems  advisable 

•  Applying  the  best  processor  configurations  obtained 
in  this  study  to  a  number  of  other  overlapping  events 
would  provide  a  firmer  estimate  of  the  achievable 
performance  levels 


V-  1 


solonos  ssrvloss  division 


Background  noise  levels  in  the  compound  events 
used  were  very  low.  The  combined  problem  of 
extracting  a  Rayleigh-  wave  signal  from  a  time¬ 
overlapping  Rayleigh-  wave  signal  and  from 
ambient  noise  needs  study 

The  possibility  of  combining  multichannel  filtering 
and  subsequent  matched  filtering  of  the  MCF  output 
may  produce  further  gains 
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the  beamsteer  processor  in  suppressing  the  interfering  event.  Multichannel  fil¬ 
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is  the  most  effective  processing  method  of  those  considered.  The  nonconvention¬ 
al  adaptive  processor  is  not  effective  in  solving  this  problem. 
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